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ABSTRACT 



We present new WFPC2 narrow band imaging of the blue compact dwarf 
galaxy I Zw 18, which is host to the lowest-metallicity H II regions known. Im- 
ages at Ha and H/3 are combined with archival broad band images to allow the 
study of the ionized gas distribution and morphology. Analysis of the Ha/H/3 
flux ratio reveals significant enhancements in some areas of both the "Northwest" 
and "Southeast" regions of the galaxy, with ratios elevated to levels as high as 
3.4. The Ha/H/5 ratio varies considerably with position throughout the galaxy. 
Comparing this distribution with the stellar distribution indicates that the re- 
gions of enhanced Ha/H/5 ratio are not due to the effects of either collisional 
excitation or underlying stellar absorption, and therefore are most likely inter- 
preted as the presence of dust. This dust has an estimated mass of (2-5) xlO 3 
M Q , which is consistent with the IRAS far-IR non-detection. Under the assump- 
tion that dust traces the presence of molecular gas, these results suggest that the 
molecular component of the ISM of I Zw 18, which is needed to fuel its active star 
formation, is also very clumpy. Such a distribution would be consistent with the 
recent FUSE non-detections of diffuse H 2 . 
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Subject headings: galaxies: evolution — galaxies: irregular — galaxies: individ- 
ual (IZwl8) — ISM: dust, extinction 



1. Introduction 

Since its discovery (Zwicky 1966), the blue compact dwarf (BCD) galaxy IZwl8 has 
been an exemplary low metallicity star forming galaxy and has undergone remarkable ob- 
servational scrutiny. The galaxy hosts the lowest metallicity H II regions known in the 
local universe, with oxygen abundances (O/H) ~ 0.02 (O/H) (Skillman & Kennicutt 1993, 
hereafter SK93). As such, the study of IZwl8 is essential in various astrophysical are- 
nas, including dwarf galaxy evolution, primordial helium abundance measurements, and star 
formation mechanisms. 

With an H I heliocentric radial velocity of 749 km distance of ~ 10 Mpc is usually 

assumed for IZwl8 (e.g., van Zee et al. 1998). Ostlin (2000) suggests a slightly larger value 
of ~ 12.6 Mpc to account for Virgocentric infall. The precise distance determination of 
I Zw 18 is beyond the scope of this paper; we follow Ostlin (2000) in adopting a distance of 
12.6 Mpc, but note that this condition is uncertain and is worthy of extensive future study. 

The low metallicity of the I Zw 18 system has led to the suggestion that it is undergoing 
its first burst of star formation (Kunth & Sargent 1986). However, the uniformity of the 
nebular abundances derived by SK93 in the two main components of the system argues 
against this hypothesis. The evolutionary status of I Zw 18 is of paramount importance, 
since the system offers the closest known reproduction of pristine conditions and therefore 
the best low metallicity environment in which to study big bang nucleosynthesis at low 
redshift. Additional physical mechanisms have also been studied in detail in IZwl8 to 
discern their behavior at very low metallicities (e.g., heating and cooling balance, Stasihska 
& Schaerer 1999). 

The nature of previous generations of star formation within I Zw 18 has not yet been 
fully characterized. Interestingly, the nebular C/O ratio was found to be elevated in IZwl8 
with respect to other comparably metal-poor irregular galaxies and above the predictions of 
massive star nucleosynthesis (Garnett et al. 1997, but see also Izotov & Thuan 1999). The 
implication of such an observation is that IZw 18 may have undergone a previous episode of 



1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc. under 
NASA contract No. NAS5-26555. 
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star formation; this conclusion is supported by the discovery of a population of 0.5-5 Gyr old 
Asymptotic Giant Branch (AGB) stars from HST/NICMOS imaging (Ostlin 2000). CMD 
analysis by Aloisi, Tosi, & Greggio (1999) also suggests the presence of evolved intermediate 
mass stars in the blue loop phase which may be as low as 3-4 M , implying ages of ~ 0.3-0.6 
Gyr. Photometric investigations thus suggest that I Zw 18 has an (at least) intermediate age 
stellar population and is thus not evolutionarily young. 

IZw 18 has been the target of deep HST studies since the launch of the spacecraft, and 
numerous archival studies have also been executed. The present new images complement 
and, where possible, improve on, previous works. Imaging studies of the nebular emission 
(Ha, He il) were carried out by Hunter & Thronson (1995, hereafter HT95) and DeMello 
et al. (1998). The former analyzed the stellar distribution and derived properties of the 
ionized gas distribution (H II region characteristics, age arguments). The latter used He II 
emission to characterize the population of Wolf-Rayet stars within the galaxy. This study 
adds to these works, in that we have deeper, higher resolution Ha imaging, and combine 
this with H/9 imaging to derive reddening characteristics within the galaxy. This is the first 
imaging diagnostic of the internal extinction using HST data. Previous evidence for dust in 
I Zw 18 has come from several ground-based spectroscopic investigations wherein reddening 
is detected (SK93; Izotov & Thuan 1998; Vilchez & Iglesias-Paramo 1998; Legrand et al. 
2000). The high spatial resolution afforded by these WFPC2 images allows us to explore the 
nature and spatial distribution of the dust and to address the possibility that the elevated 
Ha/H/5 ratios arise from collisional excitation or underlying stellar absorption. 

In this study, we seek to derive the dust properties of IZw 18. The presence of dust at 
very low metallicities has implications both for the evolution of irregular and spiral galaxies, 
and also for cosmological models of structure formation in the present epoch. It is thus 
imperative that we understand all components of such systems. This paper is organized as 
follows. In § 2 we discuss the observations and their analysis. Presentation of the new Ha 
and H/3 images, and their analysis and errors, are the topic of § 3. In § 4 we discuss the 
locations of areas rich in dust, the properties of that dust distribution, and the implications 
for the molecular gas content of the galaxy which this detection affords. We summarize our 
conclusions in § 5. 



2. Observations and Data Reduction 

Images were obtained with the Wide Field/Planetary Camera 2 (WFPC2) aboard the 
Hubble Space Telescope (HST) on 1998 April 4 and 1999 February 16. In addition, archival 
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WFPC2 data from programs GO-5309 and GO-5434 were retrieved from the STScI archive 2 
to create a comprehensive dataset covering a large wavelength regime. Here, we discuss the 
Ha and H(3 images, and their relevant continua from broad band images. A summary of 
the data used is presented in Table 1. Note that both Ha images were obtained using the 
F658N filter, because the redshift of the emission line gas in IZwl8 (Ha heliocentric radial 
velocity range of 730-780 km sec" 1 , or z~0.0025; Dufour, Esteban, & Castaneda 1996) places 
the 6563 A Ha line in the range of peak sensitivity of the F658N filter. All reductions and 
analysis were performed in the IRAF/STSDAS environment 3 . 



2.1. Calibration and Registration 

The data were initially processed through the standard HST pipeline. The pipeline- 
processed images were then corrected for geometric distortion, charge transfer efficiency, 
warm pixel effects, and cosmic ray contamination. The background level of each image was 
then removed manually. The archival broad band images were flux calibrated by multiplying 
by the image header keyword PHOTFLAM (erg sec -1 cm -2 A" 1 ), defined as the mean flux 
density that produces a count rate of 1 sec -1 , and then dividing by the header keyword 
EXPTIME (length of exposure). The On-The-Fly Calibration scheme, implemented by 
STScI, ensures all header keywords contain the most appropriate values. 

For the narrow band images, the PHOTFLAM keyword was replaced by a synthetic 
value obtained by applying the SYNPHOT package; here, the photometric parameters of an 
observation passband are calculated individually. Outputs which were used in the present 
analysis included the modified photometric calibration (output parameter URESP), the filter 
width, and a correction for emission features which appear away from the peak of a filter's 
sensitivity curve. Table 2 contains URESP and PHOTFLAM values for narrow band images, 
as well as the system response efficiency corrections which were applied. Finally, all images 
were rotated and shifted into the orientation of the new WF3 H/3 image. We find our 
registration to be accurate to within ~ 0.05" in the WF3 images (1 pixel = 0.0996") and to 
within ~ 0.03" in the PC images (1 pixel = 0.0455"). 



2 http://archive. stsci.edu 

3 IRAF is distributed by the by the National Optical Astronomy Observatories, which are operated by 
AURA, Inc., under cooperative agreement with the National Science Foundation. 
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2.2. Continuum Subtraction 

For narrow band photometry, the removal of the underlying stellar continuum is a crucial 
component of the data analysis. The contribution to the narrow band images from underlying 
continuum sources was found in a recursive manner. First, the narrow band image, which 
includes continuum contamination, was scaled by a factor roughly representative of the ratio 
of emission contained in the narrow bandpass to that contained in the total wide bandpass. 
Typically, this scaling factor was of the order of a few percent. The resultant scaled image 
was then subtracted from the wide band image whose filter response curve most closely 
isolated the redshifted emission features in question, and which contained the fewest number 
of other significant emission lines. Of course, this problem is unavoidable for certain of the 
emission line images, for example, [O III] and B./3. In such cases, the narrow band images 
were combined prior to the iterative subtraction process. This procedure was repeated by 
using different scaling factors until the wide band image no longer retained significant diffuse 
structure which is indicative of emission line radiation. Once this iterative procedure was 
completed, we were left with a pure continuum image, free from the line emission in question. 
This scaled continuum image was then subtracted from the original narrow band image(s) 
to produce a pure emission line (i.e., continuum-subtracted) map. In Figures la and lc, we 
show the PC Ha image and the underlying stellar continuum which was removed from it. 

The continuum scaling fraction was independently calculated using the SYNPHOT pack- 
age. The throughputs of the narrow band and wide filters were found at the wavelengths of 
the redshifted emission lines. The ratio of throughputs was then multiplied by the ratio of 
effective filter widths. The fractions found by these two methods agree to better than one 
percent for both emission line images. In order to obtain the total amount of line emission 
from IZwl8, the pure emission line images were then multiplied by the FWHM of each fil- 
ter, as calculated by SYNPHOT. We thus have an absolute calibration for the narrow band 
images. The continuum-subtracted images in Ha and H/5 are compared in Figures la and 
lb, respectively. The total Ha flux from I Zw 18 is found to be 3.26 x 1CT 13 erg sec -1 cm~ 2 , in 
excellent agreement with the value calculated by DeMello et al. (1998), 3.3xlCT 13 erg sec -1 
cm -2 . 



3. Emission Line Images 

The ionized gas distribution of I Zw 18 is well-resolved in the PC Ha image (Figure la). 
A wispy, filamentary structure with prominent shells dominates the morphology of the galaxy 
in the light of Ha. The very bright Ha concentration in the NW region is slightly offset from 
the peak of continuum emission (compare Figures la and lc). The Southeast (SE) region of 
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the galaxy shows more bright knots in Ha than does the NW region, and less filamentary 
structure. We do detect a smaller shell in the SE region, however; as will be shown in § 4, 
this area is also coincident with a prominent dust concentration. The relative shell sizes in 
the NW and SE regions and the presence of more OB associations and less extended diffuse 
emission in the SE region is in good agreement with the observations by HT95 and Dufour et 
al. (1996) that this region contains younger stars than the NW component. The bright H II 
regions in I Zw 18 as revealed in the PC image are shown in Figure 2 and detailed in Table 3 
in order of decreasing Ha luminosity. Our catalog compares well with that of HT95; all of 
our features are also detected and cataloged therein, except for SE 13, the least-luminous 
H II region detected, which was not visible in their F656N filter Ha image. 

The continuum-subtracted H/3 image (Figure lb) shows a very similar structure to the 
Ha image. The signal-to-noise here is lower than in the Ha image (due to lower system 
throughput than the F658N filter (roughly a factor of 2 less), a shorter integration time, and 
the intrinsically fainter B.(3 emission line strength), and produces the limiting uncertainty in 
the calculated reddening values across the galaxy (see below). Note the general agreement 
in morphology of the two images; every filament of emission in Ha is also detected in H[3, 
down to the limit of the resolution in the latter image. Resolution of fainter H II regions is 
lost in the H/3 image (e.g., SE 13; see Figures lb & 2, and Table 3). 

3.1. Reddening in I Zw 18 

Regions with high photon counts (i.e., large equivalent widths, hereafter denoted by 
EW(Ha) and EW(H/3) for Ha and H/5, respectively) in both Ha and H/5 were next analyzed 
using aperture photometry to determine fluxes. Apertures were placed on identical regions 
in six different images: Ha and H/3 fluxes, Ha and 11(3 photon counts, and continuum images 
used for each filter. To make sure that we were obtaining proper flux calibration of the Ha 
images, we performed this procedure for both the WF3 and PC F658N filter images, with 
continuum images selected to match the original resolution (i.e., PC continuum image used 
for PC Ha continuum subtraction, and likewise for the WF3 image). Table 2 (last two 
columns) contains data on the continuum subtraction values used to determine aperture flux 
ratios. 

Upon preliminary analysis of the data, we found that using small circular apertures 
centered on the H II regions (i.e., the apertures of the H II regions as shown in Figure 2 and 
listed in Table 3) yielded unacceptably large errors for most of the catalogued H II regions; 
the sources of these errors are discussed below. For this reason, we chose to expand the 
apertures to larger polygons which included the majority of visible flux when overlaid on 
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each image, encircled areas of only high EW(Ha) and EW(H/3) (compare Figures 3a and 
Id), and often contained multiple H II regions in a single aperture. This technique minimizes 
potential errors caused by slight mis-registration of the images, and enhances the signal-to- 
noise ratio by increasing photon counts in each aperture (see below for further discussion). 
The locations of these larger polygonal apertures are shown in Figure 3a. 

Flux ratios were then calculated by dividing Ha fluxes by H(3 fluxes in each aperture. 
In the absence of dust, the Ha/H/3 ratio is confined to a narrow range of values that is 
relatively insensitive to temperature and density. Following SK93, we adopt an electron 
temperature of 19,000 K for IZwl8; we thus expect Ha/H/3 = 2.76, based on the Balmer 
line ratios from Hummer & Storey (1987). Values that are elevated above this theoretical 
ratio are indicative of the presence of dust, which preferentially absorbs and scatters shorter 
wavelength radiation. We list in Table 4 the applied apertures, Ha and H/3 fluxes, and flux 
ratios with errors. We note that all Ha/H/3 ratios are consistent with the expected minimum 
theoretical value, and that the small errors make the detections of dust significant. 



3.2. Error Sources 

A potential major source of uncertainty in extragalactic reddening studies is the amount 
of foreground extinction due to gas and dust within the Galaxy. At a Galactic latitude of 
+45°, foreground extinction is expected to be low in IZwl8. In the direction of IZwl8, 
the (0?3x0?6 resolution) maps of Burstein & Heiles (1982) indicate a foreground extinction 
of E(B-V) ~ 0.01 magnitudes. The more recent, higher resolution (~ 6' FWHM) study of 
Schlegel, Finkbeiner, & Davis (1998) indicates that a marginally higher value of E(B-V) ~ 
0.03 magnitudes is applicable. Inspection of the extinction maps of this latter study indicates 
that there is no appreciable gradient in or fluctuation of E(B-V) across the luminous extent 
of I Zw 18; the average value applies over scales much larger than the field of view under 
study here (e.g., Figures l(a)-(d)), and the foreground extinction is smooth on these scales. 

For the purposes of this study, we seek to constrain the amount of internal extinction 
within IZwl8. While there is a moderate amount of extinction between the observer and 
the galaxy, which we correct for in the following analysis, any variations in the Ha/H/3 
flux ratio observed are likely to be intrinsic to IZwl8. Although we cannot rule out small- 
scale (i.e., arcsecond scales) local maxima in the foreground extinction which may affect 
our interpretation, we favor the hypothesis that searching for variations in the Ha/H/5 flux 
ratio allows us to constrain the nature of dust within I Zw 18 itself. The studies by Burstein 
& Heiles (1982) and Schlegel et al. (1998) are consistent with values of < E(B-V) < 
0.03; we adopt an intermediate value of E(B-V) = 0.02. Assuming Ry = 3.1, we then 
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have a foreground extinction of 0.06 magnitudes, which is small compared to our regions of 
highest detected extinction. Following both of the aforementioned studies, we attach to this 
foreground value a 25% uncertainty, and propagate this through our final extinction values 
calculated in Section 4.1.1. 

Uncertainties were calculated for the observed reddening values by applying Poisson 
statistics. We discuss other potential sources of error below; however, this should represent 
an accurate measure of the final error in the Ha/H/9 ratio due to the lower signal-to-noise 
ratio of the H/5 WF3 image compared to the Ha image. This dominant error term limits 
which sections of the image, and hence what physical conditions, may be sampled for the 
presence of dust. For example, the central region of the NW shell cannot be tested for the 
presence of dust, due to the fact that the number of detected H/3 photons is too small to 
allow physically meaningful errors to be assigned. 

Other sources of error in the flux calibration and image analysis procedures are discussed 
and deemed negligible as follows. We implicitly assume that removal of DC levels from images 
(e.g., background removal, scaling of images by numerical fractions) does not affect the 
standard deviation of the image, and therefore does not affect the error terms derived. Errors 
in continuum subtraction fractions were assessed by performing photometry on images with 
different continuum fractions than the final ones used. Changing the continuum subtraction 
fractions is a second order effect compared to the errors in number of counts. The values 
which we apply for continuum subtraction are also supported by the fact that the derived 
Ha/H/3 flux ratios are all consistent with the lower limit theoretical ratio of 2.76 (Hummer & 
Storey 1987, and arguments below). Furthermore, inspection of the continuum images (after 
removal of the narrow band, ionized component) reveals very little ionized gas remaining, 
indicating successful continuum subtraction (see Figures la and lc). 

We assume that the position angle of the WFPC2 camera is exactly known. Rotation 
of the images into the orientation of the new WF3 images was accomplished using simple 
linear interpolation between pixels. This method was favored due to its simplicity; we found 
no reason to apply more complicated interpolation schemes. Shifting the rotated images 
to align with the narrow band images is not an error-free process; as stated previously, we 
expect that image registration is accurate only to within a fraction of a pixel. Using the 
aforementioned larger polygonal apertures (see Figure 3a), we obtain higher photon counting 
statistics, ensuring that the contribution to the final error terms due to image rotation and 
mis-registration are of second order or higher in all cases. 
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4. Dust In I Zw 18 

In § 3, we unambiguously detect elevated Ha/H/3 flux ratios in various areas inside 
IZwl8 (see Table 4). In this section, we characterize the dust content and denote some 
implications for the evolutionary status and molecular gas content of IZwl8. 

4.1. Dust Properties 

4-1.1. Extinction, Errors And Potential Systematic Errors 
We calculate reddening values for each aperture discussed in § 3 by applying the relation 

Ih/3 Ih/3 

where I\/Ih/3 is the observed Ha/H/3 flux ratio, I\ q /Ih/3 is the theoretical ratio (2.76, 
Hummer & Storey 1987) for the adopted temperature of 19,000 K (SK93), c is the logarithmic 
reddening correction at H/5, and /(A) — f(Hj3) (= -0.37) is calculated from the interstellar 
extinction curve of Seaton (1979) as specified by Howarth (1983). We then find c and solve 
for the extinction in magnitudes by applying the conversion factor between c and Ay of 
Schild (1977), where A v = 2.17 • c. Finally, we correct for 0.06 magnitudes of foreground 
extinction (Burstein & Heiles 1982; Schlegel et al. 1998) as discussed in § 3.2. These values 
are also listed in Table 4; note that all extinctions are rather small, the largest being ~ 
0.5 magnitudes, and a fiducial value being 0.2 magnitudes for the SE region. Thus, large 
spatial variations in internal extinction in I Zw 18 are not sufficient to compromise previous 
nebular abundance studies; however, the statistically significant elevation of some Ha/H/3 
flux ratios indicates that appreciable amounts of dust do exist within the galaxy and merit 
further investigation. 

A potential difficulty in interpreting elevated Ha /H/3 ratios as indicative of the presence 
of dust is the contribution from collisional excitation of H I. Davidson & Kinman (1985) 
and SK93 have discussed the possibility that the high electron temperatures in IZw 18 may 
lead to collisional excitation of the neutral hydrogen, which would affect the emission line 
ratios of some of the lower-level hydrogen transitions, most notably Ha and B./3. Since this 
effect is directly proportional to the neutral hydrogen fraction, an observationally challenging 
quantity to measure, the importance of this effect has been difficult to assess. Photoionization 
models by Stasihska & Schaerer (1999) suggest that collisional excitation effects may affect 
the Ha/H/3 ratio by as much as 11% in the NW region of IZw 18. The distribution of 
enhanced Ha/H/9 ratios observed here argues against the importance of collisional excitation 
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of the neutral hydrogen. Most of the elevated Ha/H/3 ratios are seen in the SE component, 
which has a lower electron temperature than the NW region (SK93; Izotov & Thuan 1998). 
Furthermore, there is no evidence of enhanced Ha/H/3 flux ratios in areas of expected lower 
ionization parameter (i.e., radial gradients around ionizing stars). Thus, we interpret the 
regions which demonstrate elevated Ha/H/3 ratios (see further discussion below) to be regions 
with enhanced extinction. While we cannot rule out the possibility that collisional processes 
are making a limited contribution to the Ha emission in areas near UV-bright clusters (e.g., 
sections in the SE region; see Figure 3a), the effect is minimal even in these regions. For 
example, apertures SE D3 and SE D4 (see Table 4 and Figure 3a) are spatially separated 
by only 21.7 pc. The statistically different Ha/H/3 flux ratio between the two suggests that, 
if collisional excitation is taking place, its effects are not severe enough to compromise the 
interpretation of the photometry of this study. 

A second effect which may elevate Ha/H/3 ratios above the theoretical value is an error 
in the assumed electron temperature, T e . If T e is much lower than the adopted value, any 
elevated Ha/H/3 ratios will be caused by the preferential excitation of the lower-energy Ha 
transition compared to H/3. SK93 find T e to be lower in the SE region than in the NW 
(17,200 ± 1,200 K and 19,600 ± 600 K, respectively). However, as will be shown below, we 
find Ha/H/3 flux ratios (corrected for foreground extinction) which are as high as 3.4. To 
attain such ratios as a consequence of an overestimated T e , the actual temperature would 
have to be as low as 2,500 K in some areas, and would have to vary on unphysically short 
spatial scales (of order parsecs). Neither of these effects are observed, nor are they predicted 
by models of the photoionization structure of the galaxy (Stasihska & Schaerer 1999). We 
thus conclude that errors in the adopted T e provide a negligible error contribution in the 
following analysis. 

Finally, it is possible that underlying stellar absorption may affect the pure emission line 
images. In regions with low EW(Ha) and EW(H/3), underlying stellar Balmer absorption 
may lower the Hf3 emission line flux relative to Ha, simulating the effect of dust. To look for 
this effect we have tabulated the EW(H/3) and EW(Ha) for each of our measured regions; 
Figure Id shows the smoothed EW(Ha) map. If underlying absorption were important, 
one would expect an anti-correlation between the measured extinction and the EW(H/3). 
In fact, the data in Table 4 do show weak evidence of this trend. A regression of Ay on 
EW(H/3) yields a correlation coefficient of -0.5. This is not direct evidence that some of the 
enhancement in the Ha/H/3 ratio is due to underlying stellar absorption (it may be that 
dust is found preferentially associated with lower surface brightness stellar populations), but 
the effect does need to be considered. From the modeling of Gonzalez Delgado, Leitherer, 
& Heckman (1999), the EW(H/3) of an instantaneous burst rises monotonically from ~ 3 
A to ~ 10 A for ages of to 200 Myr. In a constant star formation model EW(H/3) rises 
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monotonically from ~ 3 A to ~ 6 A over the same period. If we could be confident that all 
of the stars in IZwl8 were younger than 20 Myr, then an upper limit of 6 A EW(H/3) of 
underlying H/3 absorption would be reasonable. If there are stars older than this in IZwl8, 
then the continuous star formation model would be more appropriate (which has a maximum 
value of ~ 8 A EW(H/9) for ages of ~ 0.5 Gyr and higher). For our lowest EW(H/9) regions 
(~75 to ~100 A), some, if not most, of the enhancement in the Ha/H/3 ratio could be due 
to underlying stellar absorption if it were mixed with an intermediate age stellar population. 
However, for the regions with A v in excess of 0.15 (Ha/H/3 = 3.00) it is unlikely that the 
enhancement in the Ha/H/3 ratio could be due solely to underlying stellar absorption. Note 
that in several regions (with EW(H/3) = 120 to 160 A, see Table 4) the Ha/H/3 ratio is 
consistent with no reddening (and thus no effect of underlying absorption). Thus, while we 
cannot completely rule out some enhancement in the Hct/H/3 ratio due to underlying stellar 
absorption, we can rule out that it is the dominant effect in the regions where we have the 
highest Ha/H/5 ratios. Since these regions are the focus of what follows, we will assume that 
the effects of underlying absorption are negligible. 



4-1.2. Dust Locations And Masses 

In Figure 3b, we overlay apertures on the PC Ha image of IZwl8 which outline the 
locations of elevated Ha/H/3 ratios (note that these trapezoidal apertures are meant to guide 
the eye to areas with Ha/H/9 > 3.0, and in some cases contain more than one of the apertures 
used to derive the Ha/H/3 flux ratios). We detect dust in regions which are near UV-bright 
sources, as well as in regions which are far away from them. The high surface brightness H II 
regions which have significant extinction are all located in the SE component, in agreement 
with the ground based optical observations (SK93; Izotov & Thuan 1998; Vilchez & Iglesias- 
Paramo 1998). This is expected, since the SE region contains younger stars and likely 
more dense molecular clouds than the NW region. The distribution of dust (and therefore 
molecular gas) suggests that future star formation is likely in the SE region; furthermore, the 
reddening morphology hints at the presence of a band of extinction between the SE and NW 
components. This region cannot be sampled for dust with the present Ha/H/9 ratio method, 
due to the absence of Ha and H/3 emission in this region (see § 3.2). Overall, we detect dust 
with a similar distribution to that found in the dwarf starburst galaxy NGC 5253 (Calzetti 
et al. 1997). The irregular nature of the dust distribution immediately implies one or both 
of the following scenarios. First, dust which is near UV-bright clusters is surviving in a very 
hard radiation field, implying an efficient mechanism of self-shielding on small scales and the 
prominence of clumpy structure in the ISM. Second, dust which is displaced from UV-bright 
clusters is consistent with dust production by a previous generation of star formation. 
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In order to estimate a mass of the dust which we detect, we must make certain as- 
sumptions about the nature of the grains. One method with which to obtain a dust mass 
estimate is to use the ratio of intensities of far infrared (FIR) emission at 60 and 100 /zm. 
However, I Zw 18 is not detected by the IRAS satellite in either waveband. A second method 
with which to find the dust mass is to apply a detailed balance argument, assuming one 
knows the nature of the incident radiation field and the composition of the grains. How- 
ever, especially in the case of BCD's such as IZwl8 where the nature of the underlying 
stellar population is not yet fully characterized, this detailed balance argument is also not 
applicable. 

A third and widely used method to obtain a dust mass is to apply a gas-to-dust ratio 
which is applicable to similar physical conditions (metallicity, radiation field, star formation 
efficiency, etc.). Seminal Galactic work on this correlation by Bohlin, Savage, & Drake (1978) 
finds that the mean ratio <(H I + H 2 )/E(B-V)> is nearly constant (within a factor of 1.5), 
with a value of 5.8 xlO 21 cm -2 mag" 1 . From a larger stellar sample, Diplas & Savage (1994) 
find a slightly lower value of 4.93 xlO 21 cm" 2 mag" 1 . Extinction studies in other environ- 
ments suggest a trend of increasing gas-to-dust ratio with decreasing ambient metallicity. 
Extinction studies toward stars in the LMC arrive at an average value of <(Nh i)/E(B-V)> 
= 2.4 xlO 22 cm" 2 mag" 1 (Fitzpatrick 1985; Dufour, Shields, & Talbot 1982; Clayton et al. 
1996). The more metal-poor conditions of the SMC yield yet higher values, perhaps as high 
as 10 23 cm" 2 mag" 1 (Prevot et al. 1984). 

The H I synthesis imaging study of van Zee et al. (1998) shows that there is a column 
density maximum of 3xl0 21 cm" 2 which is coincident with the SE component (see their 
Figure 6). If we assume that the SE ionized component is embedded in an H I envelope, 
then we can make a very crude estimate of the ratio of H I column density to extinction by 
applying the relation 

N(HI) 1.5 xlO 21 cm" 2 n n ~ 2 , 

~ 2.3 x 10 22 cm" 2 mag" 1 (2) 



E(B - V) 0.06 mag 

where we have assumed Ry = Ay/E(B-V) = 3.1, and that a representative value of the 
average extinction towards the SE region is Ay = 0.2 mag. This value is comparable to, 
but less than, estimates of the gas-to-dust ratio in the SMC (3.7-5.2xl0 22 cm" 2 mag" 1 ; 
Bouchet et al. 1985). However, it makes no correction for a potentially large contribution 
from molecular hydrogen (necessary to sustain the high star formation rate; see arguments 
below), which may allow the two results to be consistent. Furthermore, this estimate only 
accounts for the fraction of the dust which may be foreground to I Zw 18; if a substantial 
amount of dust is located behind or mixed within the H II regions, this argument will (perhaps 
significantly) underestimate the total dust mass within the galaxy. 
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Next, we must assume a value for the dust-to-gas ratio by mass in order to obtain a dust 
mass estimate. We make such an assumption in two ways. First, we note that Stanimirovic 
et al. (2000) finds the average dust-to-gas ratio in the SMC to be 8.2xl0~ 5 by mass (roughly 
30 times less than in the Galaxy). Since this is the lowest-metallicity galaxy which has 
well studied extinction properties, we naively apply this ratio to IZw 18. This represents an 
apparent upper limit, since there is evidence of a trend towards lower dust-to-gas ratios at 
lower metallicities, as previously mentioned, and the metallicity of I Zw 18 is ~ 0.8 dex lower 
than that of the SMC. Applying such an assumption and Equation 2, rewritten in the more 
convenient form 

A V /N H = 1.6 x 10~ 22 mag cm 2 (3) 

we arrive at a dust mass of ~ 2000 M© (with no correction for the molecular gas content). 

The second method with which to estimate the dust-to-gas ratio is to apply the analytical 
relation found by Lisenfeld & Ferrara (1998, hereafter LF98), 12 + log(0/H) oc (0.52±0.25) • 
log(M dust /M H i). While this relation was derived for dwarf irregulars and not BCD's, we 
include it for completeness, and because the application of the relation to the SMC yields a 
dust-to-gas ratio of 6.5 xlO" 5 by mass, in general agreement with the observational result of 
Stanimirovic et al. (2000). Applying the same relation to IZw 18, we obtain a dust-to-gas 
ratio of 1.7xl0~ 4 by mass, which is consistent with the upper limit of 3.7x10 obtained by 
LF98 by assuming a 100 /im flux for I Zwl8. We thus assume a value of 2xl0~ 4 for the dust- 
to-gas ratio using this relation, but note its general uncertainty, and that it appears high in 
comparison to observational estimates from other low- metallicity systems (e.g., Stanimirovic 
et al. 2000). Application of Equation 3 then yields a dust mass of ~ 4800 M . 

From the above arguments, we estimate a total dust mass of (2-5) xlO 3 M within 
IZw 18, due to the lack of further constraints which may be placed on the nature of the 
dust content from these observations. We note that this dust mass is lower than all of the 
BCD's which are detected in the IRAS bands in the study by LF98, as expected for the very 
low metallicity of IZw 18. We thus consider this a conservative estimate of the total dust 
content of this BCD. If the dust production process and evolution behaves as in the dwarf 
irregular sample (e.g., all dwarf irregulars with 12 + log(0/H) < 7.4 have dust masses < 
10 3 M Q ; LF98), then the above values may be too high. Further constraints on the nature 
of the dust content of I Zw 18 will be very useful. 

4.2. Implications For The Molecular Gas Content 

I Zw 18 is undergoing an intense star formation episode, as evidenced by its blue photo- 
metric color (U— B = -0.88, B— V = -0.03; van Zee et al. 1998), its large number of ionizing, 
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massive stars (see Table 3), and by the large ionized gas extent (e.g., Figures la and lb). 
Consequently, we expect concentrations of molecular gas throughout the galaxy, since we 
observe large H II regions and OB associations. I Zw 18 is not detected in CO observations, 
nor is any galaxy with ambient metallicity below 12 + log(0/H) ~ 8.0 (Taylor, Kobul- 
nicky & Skillman 1998; Barone et al. 2000). The potential reasons why detections of CO 
in low-metallicity galaxies have remained elusive are numerous. One explanation is that 
the conversion factor, X co , from I co to N Ha becomes nonlinear below a certain metallicity 
(Maloney & Black 1988). At low abundances, the column density of the CO molecule will 
decrease. In areas undergoing massive star formation, the hard UV radiation field will more 
easily photodissociate CO to the point where it can no longer be self-shielding; Madden et 
al. (1997) find observational evidence for such behavior. This will decrease the component 
of molecular clouds which will have prominent CO emission, while the H2 is less affected. 
The filling factor decreases, reducing the CO luminosity for a given molecular gas mass, and 
incompletely sampling the H 2 mass as a result. 

These large dust concentrations (of order 50-100 pc, see Figure 3b) suggest that the 
process of H2 formation on grain surfaces may be ongoing in IZwl8. Recent FUSE non- 
detections of H 2 in absorption in IZwl8 (Vidal-Madjar et al. 2000) can still be consistent 
with the presence of large amounts of molecular gas, since these HST results indicate that at 
least some of the dust is peripheral to the UV-bright clusters, and not necessarily along lines 
of sight to them. We interpret these results as evidence for a clumpy ISM in I Zw 18; small, 
dense molecular clouds will allow active star formation, while their volume filling factor and 
overall CO luminosity will be very small. Vidal-Madjar et al. (2000) also conclude that 
dense, discrete clouds may have gone undetected in their FUSE observations. Clearly, very 
high-resolution observations (i.e., interferometric CO mapping; e.g., Walter et al. 2001) of 
the ISM of I Zw 18 and of other low-metallicity galaxies are needed to attempt to locate the 
molecular gas. 

It is most intriguing that the two lowest-metallicity galaxies known, I Zw 18 and SBS 0335- 
052 (12+log(0/H) = 7.33±0.01, Izotov et al. 1997), both contain significant amounts of dust. 
As discussed above, in I Zw 18, some of the dust concentrations are displaced from ionized 
regions. Similarly, in SBS 0335-052, dust is found displaced from the locations of the six 
massive star clusters (Thuan, Izotov, & Lipovetsky 1997). In both galaxies, the presence of 
dust displaced from some of the ionizing clusters suggests that the molecular gas content of 
low-metallicity galaxies is very difficult to quantify via UV absorption or CO tracer lines. In 
particular, these results suggest that the lack of molecular gas as inferred from such investi- 
gations does not necessarily imply a lack of molecular gas, but rather a different mechanism 
of shielding (e.g., small, dense molecular clouds) when the metallicity is very low. 
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5. Conclusion 

New WFPC2 narrow band images of the blue compact dwarf galaxy IZw 18 have been 
analyzed. Images at Ha and Hf3 are used to constrain the dust content and distribution of 
the galaxy to analyze the morphology of the ionized gas, and to characterize individual H II 
regions. We detect significantly elevated Ha/H/3 flux ratios throughout the galaxy, both close 
to and displaced from UV sources, and in general agreement with the morphology derived by 
Vilchez & Iglesias-Paramo (1998). The effects of collisional excitation and stellar absorption 
are shown to be minimal, and we interpret these ratios as indicative of the presence of dust. 
The small total dust mass, estimated to be of order 10 3 M Q , is expected for the very low 
metallicity of the galaxy and is consistent with the IRAS non-detection in the far infrared 
regime. The presence of dust at low metallicities is now seen in the two lowest metallicity 
galaxies and is significant, since these may be the most ideal testbeds for studies of galaxy 
formation and chemical evolution processes at low redshifts. 

The presence of dust in I Zw 18 allows us to form a consistent picture of the galaxy's 
interstellar medium as derived from various observational studies. I Zw 18 is a low mass 
galaxy which is embedded in a large H I envelope (van Zee et al. 1998). Within this neutral 
gas, we see in deep optical and infrared imaging (HT95; DeMello et al. 1998; Aloisi et al. 
1999; Ostlin 2000) the evolutionary signatures of stars down to the lower mass limit which 
is currently detectable (~ 1.1 M Q in NICMOS imaging; Ostlin 2000). The large ionized 
regions and intense Ha emission evidence the UV flux of massive stars and suggest ongoing 
star formation. Mixed throughout the galaxy we find large dust concentrations, which could 
be due to either (or both) mass loss from massive stars or intermediate mass (AGB) stars that 
have been recently supported by the aforementioned CMD studies. Regardless of its origin, 
however, the presence of dust suggests that areas rich in molecular gas are abundant even at 
this very low metallicity. The FUSE non-detection of diffuse H 2 in absorption (Vidal-Madjar 
et al. 2000) may still be consistent with the presence of large amounts of molecular gas in 
I Zw 18, since in these observations we find bulk dust distributions displaced from UV clusters 
and not along lines of sight directly to them. We favor the scenario that the molecular clouds 
in I Zw 18 (and likely in other low-metallicity galaxies) are small and dense. This will cause 
them to be very difficult to detect, yet still allow them to actively and efficiently form the 
short-lived massive stars which we see throughout I Zw 18. 
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Fig. I. — Continuum-subtracted Ha (a) and H/3 (b) pure emission line images of IZwl8, 
compared with the scaled emission line subtracted R-band image (c) and the (smoothed) 
equivalent width of the Ha emission (d). The scale encompasses the full range of intensity 
for all images, except for (d), where the scale is from 100 A or less (white) to 2000 A or 
greater (black). The field of view, identical in all images, is 13.7" x 10.5", or 833 pc x 641 
pc at the adopted distance of IZwl8. The arrows denote north (tipped) and east and are 
1" in length. The lower signal-to-noise ratio and resolution of the H/3 image compared to 
the Ha image will form the dominant contribution to the errors in the Ha/H/3 flux ratio. 
Note the negligible residual ionized gas in the stellar continuum, which is slightly displaced 
from the ionized gas; there are also stars which are outside the luminous extent of IZwl8 
in Ha. The Ha equivalent width map (d) was smoothed to ~ 0.25" resolution, and was 
masked to only retain information in areas of high signal-to-noise ratio in the original Ha 
image (for clarity). A comparison of (d) with the aperture placements used to derive the 
Ha/H/? flux ratios (see Figure 3a, and discussion in Sections 3.1 and 4.1.1) shows that all 
apertures encompass areas with large equivalent width (> 700 A in Ha and > 80 A in H/5; 
see Table 4), making negligible the potential effects of stellar absorption. 
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Fig. 2. — The bright H II regions of I Zw 18, overlaid on the continuum-subtracted PC Ha 
image. The characteristics of the regions, numbered in order of decreasing Ha luminosity, 
are summarized in Table 3. 
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Fig. 3. — Continuum-subtracted Ha image of I Zw 18, overlaid with the locations of apertures 
used to determine Ha/H/3 flux ratios (a). The apertures are labeled, and the reddening and 
extinction to each is shown in Table 4. In (b), the same image is overlaid with the locations 
of significant dust concentrations (i.e., Ha/H/5 > 3.0). Some of the dust-rich regions are near 
ionizing sources (e.g., in the SE region), while others are displaced from clear sources of UV 
flux. These dust concentrations may isolate areas which are rich in molecular gas; note that 
these are not necessarily coincident with sightlines to UV-bright clusters (c.f., Vidal-Madjar 
et al. 2000). 
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Table I. HST WFPC2 Observations of I Zw 18 



Filter 


Emission 


Proposal 


Date 


Exposure 


WFPC2 


Dataset 










(sec) 


Chip 




F487N 


Hp 


GO-6536 


2/16/99 


2500 


WF3 


U39A0101R,U39A0102R,U39A0103R 


F658N 


Ha 


GO-6536 


4/04/98 


2500 


PC 


U39A0201R,U39A0202R,U39A0203R 


F658N 


Ha 


GO-5434 


11/02/94 


4600 


WF3 


U2F90205T,U2F90206T 


F450W 


Wide-B 


GO-5434 


11/03/94 


4600 


WF3 


U2F90102T,U2F90103T 


F675W 


R 


GO-5434 


3/01/95 


2000 


PC 


U2F90305T,U2F90306T 


F702W 


Wide-R 


GO-5434 


11/03/94 


5400 


WF3 


U2F90101T,U2F90201T,U2F90202T 
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Table 2. Narrow Band Photometric Calibration Values 



Filter 


Emission 


WFPC2 


PHOTFLAM a 


URESP a 


Filter 


Efficiency 


Continuum 


Scaled b 






Chip 






Width (A) 


Correction 


Filter 


Fraction 


F487N 


H/3 


WF3 


3.89xl0~ 16 


3.94xl0~ 16 


33.9 


1.07 


F450W 


2.7% 


F658N 


Ha 


PC 


1.06xl0~ 16 


1.04X10" 16 


39.2 


1.01 


F702W 


1.9% 


F658N 


Ha 


WF3 


1.06X10" 16 


1.04X10" 16 


39.2 


1.01 


F675W 


2.9% 



a Units of erg sec 1 cm 2 A 1 

b Fraction which narrow band image was multiplied by before subtraction from wide band filter. This value was found 
both iteratively and by finding the ratio of (filter throughput at the wavelength of redshifted emission) to (filter width at 
full width half maximum); both ratios were calculated using the SYNPHOT package. Detailed discussion of the continuum 
subtraction process is given in § 2.2. 
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Table 3. The Bright H II Regions of I Zw 18 



Feature 


Center RA 


Center DEC 


Previous 


Radius 




08 Ionizing 




(J2000) 


(J2000) 


Idcntification a 


(PC) 


(10 36 erg sec" 1 ) 


Equivalent 15 


NW 1 


9:34:02.361 


55:14:26.27 


HT 2 


11.2 


511 


238 


SE 1 


9:34:02.610 


55:14:21.45 


HT 16 


7.5 


112.8 


53 


NW 2 


9:34:02.428 


55:14:26.11 


HT 3 


6.0 


71.1 


33 


SE 2 


9:34:02.680 


55:14:21.75 


HT 5 


5.0 


42.3 


21 


SE 3 


9:34:02.660 


55:14:20.35 


HT 13 


3.4 


37.7 


18 


SE 4 


9:34:02.473 


55:14:20.73 


HT 10 


4.5 


35.0 


16 


SE 5 


9:34:02.549 


55:14:22.39 


HT 6 


2.5 


29.7 


14 


SE 6 


9:34:02.680 


55:14:20.62 


HT 14 


3.0 


29.1 


14 


SE 7 


9:34:02.526 


55:14:22.13 


HT 7 


3.00 


23.0 


11 


SE 8 


9:34:02.604 


55:14:20.84 


HT 12 


3.0 


21.7 


10 


SE 9 


9:34:02.565 


55:14:20.92 


HT 11 


2.9 


18.1 


8 


SE 10 


9:34:02.452 


55:14:21.34 


HT 9 


2.4 


17.1 


8 


SE 11 


9:34:02.416 


55:14:21.49 


HT 8 


2.6 


16.2 


8 


SE 12 


9:34:02.525 


55:14:20.98 


HT 11 


2.5 


12.5 


6 


SE 13 


9:34:02.507 


55:14:21.78 




2.85 


9.2 


4 


Total Galaxy 


9:34:02.0 


55:14:24.0 






6200 


2886 



a HT = Hunter & Thronson (1995) = HT 95 

b Assumes an ionizing photon rate of 10 48 ' 6 sec -1 for an 08 star (Osterbrock 1989). 
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Table 4. Polygonal Aperture Flux Ratios and Extinction 



Feature 


Ha Flux 

(erg sec -1 cm -2 ) 


H/3 Flux 
(erg sec -1 cm~ 2 ) 


Ha/H/3 
Ratio 


A v a 
(Magnitudes) 


Ha Equivalent 
Width b (A) 


H/3 Equivalent 
Width b (A) 


NW Dl 


3.88xl(r 14 


1.34xl0~ 14 


2.88±0.02 


0.05±0.01 


908±3 


88±1 


NW D2 


7.13X10- 15 


2.49xl0~ 15 


2.87±0.04 


0.04±0.02 


2055±20 


192±3 


NW D3 


2.05xlCr 14 


6.45x10" 15 


3.17±0.03 


0.29±0.02 


903±4 


75±1 


NW D4 


3.02xlCr 14 


1.07xl(r 14 


2.83±0.02 


0.00±0.01 


1527±7 


156±1 


NW D5 


4.87X10" 14 


1.65xl0~ 14 


2.96±0.02 


0.12±0.01 


1342±4 


140±1 


NW D6 


4.17X10" 15 


1.37xl0- 15 


3.04±0.06 


0.19±0.03 


1412±15 


137±3 


SE Dl 


9.59xlCr 15 


3.20xl0~ 15 


3.00±0.04 


0.15±0.02 


920±6 


95±1 


SE D2 


5.40xl0~ 15 


1.93X1CT 15 


2.80±0.04 


0.00±0.02 


1067± 9 


126±2 


SE D3 


6.81 xlO~ 15 


2.48xl0~ 15 


2.75±0.04 


0.00±0.03 


1175±12 


150±3 


SE D4 


1.82X10" 18 


5.35xl0- lfi 


3.40±0.11 


0.47±0.05 


1255±20 


101±3 


SE D5 


9.66xl0~ 15 


3.07xl0~ 15 


3.15±0.04 


0.28±0.02 


1369±10 


115±2 


SE D6 


1.12X10" 14 


3.73xl0~ 15 


3.00±0.04 


0.15±0.02 


1110±7 


101±1 


SE D7 


5.51 xlO~ 15 


1.87X1CT 15 


2.94±0.05 


0.10±0.03 


1103±9 


103±2 


SE D8 


4.72xl0~ 15 


1.56xl0~ 15 


3.03±0.06 


0.18±0.03 


767±6 


84±2 



a Corrected for 0.06 magnitudes of foreground extinction. 
b Average value over the aperture, rounded to nearest A. 



